Integrating terrestrial and aquatic ecosystems to constrain estimates of land-atmosphere carbon exchange.

In this Perspective, we put forward an integrative framework to improve estimates of land-atmosphere carbon exchange based on the accumulation of carbon in the landscape as constrained by its lateral export through rivers. The framework uses the watershed as the fundamental spatial unit and integrates all terrestrial and aquatic ecosystems as well as their hydrologic carbon exchanges. Application of the framework should help bridge the existing gap between land and atmosphere-based approaches and offers a platform to increase communication and synergy among the terrestrial, aquatic, and atmospheric research communities that is paramount to advance landscape carbon budget assessments.

N and P constrain C in ecosystems under climate change: Role of nutrient redistribution, accumulation, and stoichiometry.

We use the Multiple Element Limitation (MEL) model to examine responses of 12 ecosystems to elevated carbon dioxide (CO2 ), warming, and 20% decreases or increases in precipitation. Ecosystems respond synergistically to elevated CO2 , warming, and decreased precipitation combined because higher water-use efficiency with elevated CO2 and higher fertility with warming compensate for responses to drought. Response to elevated CO2 , warming, and increased precipitation combined is additive. We analyze changes in ecosystem carbon (C) based on four nitrogen (N) and four phosphorus (P) attribution factors: (1) changes in total ecosystem N and P, (2) changes in N and P distribution between vegetation and soil, (3) changes in vegetation C:N and C:P ratios, and (4) changes in soil C:N and C:P ratios. In the combined CO2 and climate change simulations, all ecosystems gain C. The contributions of these four attribution factors to changes in ecosystem C storage varies among ecosystems because of differences in the initial distributions of N and P between vegetation and soil and the openness of the ecosystem N and P cycles. The net transfer of N and P from soil to vegetation dominates the C response of forests. For tundra and grasslands, the C gain is also associated with increased soil C:N and C:P. In ecosystems with symbiotic N fixation, C gains resulted from N accumulation. Because of differences in N versus P cycle openness and the distribution of organic matter between vegetation and soil, changes in the N and P attribution factors do not always parallel one another. Differences among ecosystems in C-nutrient interactions and the amount of woody biomass interact to shape ecosystem C sequestration under simulated global change. We suggest that future studies quantify the openness of the N and P cycles and changes in the distribution of C, N, and P among ecosystem components, which currently limit understanding of nutrient effects on C sequestration and responses to elevated CO2 and climate change.

Resilience and sensitivity of ecosystem carbon stocks to fire-regime change in Alaskan tundra.

Fire disturbance has increased in some tundra ecosystems due to anthropogenic climate change, with important ramifications for terrestrial carbon cycling. Assessment of the potential impact of fire-regime change on tundra carbon stocks requires long-term perspectives because tundra fires have been rare historically. Here we integrated the process-based Dynamic Organic Soil version of the Terrestrial Ecosystem Model with paleo-fire records to evaluate the responses of tundra carbon stocks to changes in fire return interval (FRI). Paleorecords reveal that mean FRIs of tundra ecosystems in Alaska ranged from centennial to millennial timescales (200-6000 years) during the late Quaternary, but projected FRIs by 2100 decrease to a few hundred years to several decades (70-660 years). Our simulations indicate threshold effects of changing FRIs on tundra carbon stocks. Shortening FRI from 5000 to 1000 years results in minimal carbon release (<5%) from Alaskan tundra ecosystems. Rapid carbon stock loss occurs when FRI declines below 800 years trigger sustained mobilization of ancient carbon stocks from permafrost soils. However, substantial spatial heterogeneity in the resilience/sensitivity of tundra carbon stocks to FRI change exists, largely attributable to vegetation types. We identified the carbon stocks in shrub tundra as the most vulnerable to decreasing FRI because shrub tundra stores a large share of carbon in combustible biomass and organic soils. Moreover, our results suggest that ecosystems characterized by large carbon stocks and relatively long FRIs (e.g. Brooks Foothills) may transition towards hotspots of permafrost carbon emission as a response to crossing FRI thresholds in the coming decades. These findings combined imply that fire disturbance may play an increasingly important role in future carbon balance of tundra ecosystems, but the net outcome may be strongly modulated by vegetation composition.

Assessing dynamic vegetation model parameter uncertainty across Alaskan arctic tundra plant communities.

As the Arctic region moves into uncharted territory under a warming climate, it is important to refine the terrestrial biosphere models (TBMs) that help us understand and predict change. One fundamental uncertainty in TBMs relates to model parameters, configuration variables internal to the model whose value can be estimated from data. We incorporate a version of the Terrestrial Ecosystem Model (TEM) developed for arctic ecosystems into the Predictive Ecosystem Analyzer (PEcAn) framework. PEcAn treats model parameters as probability distributions, estimates parameters based on a synthesis of available field data, and then quantifies both model sensitivity and uncertainty to a given parameter or suite of parameters. We examined how variation in 21 parameters in the equation for gross primary production influenced model sensitivity and uncertainty in terms of two carbon fluxes (net primary productivity and heterotrophic respiration) and two carbon (C) pools (vegetation C and soil C). We set up different parameterizations of TEM across a range of tundra types (tussock tundra, heath tundra, wet sedge tundra, and shrub tundra) in northern Alaska, along a latitudinal transect extending from the coastal plain near Utqiagvik to the southern foothills of the Brooks Range, to the Seward Peninsula. TEM was most sensitive to parameters related to the temperature regulation of photosynthesis. Model uncertainty was mostly due to parameters related to leaf area, temperature regulation of photosynthesis, and the stomatal responses to ambient light conditions. Our analysis also showed that sensitivity and uncertainty to a given parameter varied spatially. At some sites, model sensitivity and uncertainty tended to be connected to a wider range of parameters, underlining the importance of assessing tundra community processes across environmental gradients or geographic locations. Generally, across sites, the flux of net primary productivity (NPP) and pool of vegetation C had about equal uncertainty, while heterotrophic respiration had higher uncertainty than the pool of soil C. Our study illustrates the complexity inherent in evaluating parameter uncertainty across highly heterogeneous arctic tundra plant communities. It also provides a framework for iteratively testing how newly collected field data related to key parameters may result in more effective forecasting of Arctic change.

Mycobiont contribution to tundra plant acquisition of permafrost-derived nitrogen.

As Arctic soils warm, thawed permafrost releases nitrogen (N) that could stimulate plant productivity and thus offset soil carbon losses from tundra ecosystems. Although mycorrhizal fungi could facilitate plant access to permafrost-derived N, their exploration capacity beyond host plant root systems into deep, cold active layer soils adjacent to the permafrost table is unknown. We characterized root-associated fungi (RAF) that colonized ericoid (ERM) and ectomycorrhizal (ECM) shrub roots and occurred below the maximum rooting depth in permafrost thaw-front soil in tussock and shrub tundra communities. We explored the relationships between root and thaw front fungal composition and plant uptake of a 15 N tracer applied at the permafrost boundary. We show that ERM and ECM shrubs associate with RAF at the thaw front providing evidence for potential mycelial connectivity between roots and the permafrost boundary. Among shrubs and tundra communities, RAF connectivity to the thaw boundary was ubiquitous. The occurrence of particular RAF in both roots and thaw front soil was positively correlated with 15 N recovered in shrub biomass Taxon-specific RAF associations could be a mechanism for the vertical redistribution of deep, permafrost-derived nutrients, which may alleviate N limitation and stimulate productivity in warming tundra.

Spatiotemporal remote sensing of ecosystem change and causation across Alaska.

Contemporary climate change in Alaska has resulted in amplified rates of press and pulse disturbances that drive ecosystem change with significant consequences for socio-environmental systems. Despite the vulnerability of Arctic and boreal landscapes to change, little has been done to characterize landscape change and associated drivers across northern high-latitude ecosystems. Here we characterize the historical sensitivity of Alaska's ecosystems to environmental change and anthropogenic disturbances using expert knowledge, remote sensing data, and spatiotemporal analyses and modeling. Time-series analysis of moderate-and high-resolution imagery was used to characterize land- and water-surface dynamics across Alaska. Some 430,000 interpretations of ecological and geomorphological change were made using historical air photos and satellite imagery, and corroborate land-surface greening, browning, and wetness/moisture trend parameters derived from peak-growing season Landsat imagery acquired from 1984 to 2015. The time series of change metrics, together with climatic data and maps of landscape characteristics, were incorporated into a modeling framework for mapping and understanding of drivers of change throughout Alaska. According to our analysis, approximately 13% (~174,000 ± 8700 km2 ) of Alaska has experienced directional change in the last 32 years (±95% confidence intervals). At the ecoregions level, substantial increases in remotely sensed vegetation productivity were most pronounced in western and northern foothills of Alaska, which is explained by vegetation growth associated with increasing air temperatures. Significant browning trends were largely the result of recent wildfires in interior Alaska, but browning trends are also driven by increases in evaporative demand and surface-water gains that have predominately occurred over warming permafrost landscapes. Increased rates of photosynthetic activity are associated with stabilization and recovery processes following wildfire, timber harvesting, insect damage, thermokarst, glacial retreat, and lake infilling and drainage events. Our results fill a critical gap in the understanding of historical and potential future trajectories of change in northern high-latitude regions.

Large loss of CO2 in winter observed across the northern permafrost region.

Recent warming in the Arctic, which has been amplified during the winter1-3, greatly enhances microbial decomposition of soil organic matter and subsequent release of carbon dioxide (CO2)4. However, the amount of CO2 released in winter is highly uncertain and has not been well represented by ecosystem models or by empirically-based estimates5,6. Here we synthesize regional in situ observations of CO2 flux from arctic and boreal soils to assess current and future winter carbon losses from the northern permafrost domain. We estimate a contemporary loss of 1662 Tg C yr-1 from the permafrost region during the winter season (October through April). This loss is greater than the average growing season carbon uptake for this region estimated from process models (-1032 Tg C yr-1). Extending model predictions to warmer conditions in 2100 indicates that winter CO2 emissions will increase 17% under a moderate mitigation scenario-Representative Concentration Pathway (RCP) 4.5-and 41% under business-as-usual emissions scenario-RCP 8.5. Our results provide a new baseline for winter CO2 emissions from northern terrestrial regions and indicate that enhanced soil CO2 loss due to winter warming may offset growing season carbon uptake under future climatic conditions.

Assessing historical and projected carbon balance of Alaska: A synthesis of results and policy/management implications.

We summarize the results of a recent interagency assessment of land carbon dynamics in Alaska, in which carbon dynamics were estimated for all major terrestrial and aquatic ecosystems for the historical period (1950-2009) and a projection period (2010-2099). Between 1950 and 2009, upland and wetland (i.e., terrestrial) ecosystems of the state gained 0.4 Tg C/yr (0.1% of net primary production, NPP), resulting in a cumulative greenhouse gas radiative forcing of 1.68 × 10-3  W/m2 . The change in carbon storage is spatially variable with the region of the Northwest Boreal Landscape Conservation Cooperative (LCC) losing carbon because of fire disturbance. The combined carbon transport via various pathways through inland aquatic ecosystems of Alaska was estimated to be 41.3 Tg C/yr (17% of terrestrial NPP). During the projection period (2010-2099), carbon storage of terrestrial ecosystems of Alaska was projected to increase (22.5-70.0 Tg C/yr), primarily because of NPP increases of 10-30% associated with responses to rising atmospheric CO2 , increased nitrogen cycling, and longer growing seasons. Although carbon emissions to the atmosphere from wildfire and wetland CH4 were projected to increase for all of the climate projections, the increases in NPP more than compensated for those losses at the statewide level. Carbon dynamics of terrestrial ecosystems continue to warm the climate for four of the six future projections and cool the climate for only one of the projections. The attribution analyses we conducted indicated that the response of NPP in terrestrial ecosystems to rising atmospheric CO2 (~5% per 100 ppmv CO2 ) saturates as CO2 increases (between approximately +150 and +450 ppmv among projections). This response, along with the expectation that permafrost thaw would be much greater and release large quantities of permafrost carbon after 2100, suggests that projected carbon gains in terrestrial ecosystems of Alaska may not be sustained. From a national perspective, inclusion of all of Alaska in greenhouse gas inventory reports would ensure better accounting of the overall greenhouse gas balance of the nation and provide a foundation for considering mitigation activities in areas that are accessible enough to support substantive deployment.

Functional responses of white spruce to snowshoe hare herbivory at the treeline.

Herbivores can modify the rate of shrub and treeline advance. Both direct and indirect effects of herbivory may simultaneously interact to affect the growth rates of plants at this ecotone. We investigated the effect of snowshoe hare herbivory on the height of white spruce at two treeline locations in Alaska, USA. White spruce is expanding its distribution both upwards in elevation and northward in latitude because of climate warming, and snowshoe hares are already present in areas likely to be colonized by spruce. We hypothesized that herbivory would result in browsed individuals having reduced height, suggesting herbivory is a direct, negative effect on spruce treeline advance. We found an interactive effect between browsing history and spruce age. When young (under 30 years old), individuals that were browsed tended to be taller than unbrowsed individuals. However, older seedlings (over 30 years old) that had been browsed were shorter than unbrowsed individuals of the same age. Hares suppress faster growing individuals that are initially taller by preferentially browsing them as they emerge above the winter snowpack. This reduced height, in combination with increased mortality associated with browsing, is predicted to slow the advance of both latitudinal and altitudinal treeline expansions and alter the structure of treeline forests.

The role of environmental driving factors in historical and projected carbon dynamics of wetland ecosystems in Alaska.

Wetlands are critical terrestrial ecosystems in Alaska, covering ~177,000 km2 , an area greater than all the wetlands in the remainder of the United States. To assess the relative influence of changing climate, atmospheric carbon dioxide (CO2 ) concentration, and fire regime on carbon balance in wetland ecosystems of Alaska, a modeling framework that incorporates a fire disturbance model and two biogeochemical models was used. Spatially explicit simulations were conducted at 1-km resolution for the historical period (1950-2009) and future projection period (2010-2099). Simulations estimated that wetland ecosystems of Alaska lost 175 Tg carbon (C) in the historical period. Ecosystem C storage in 2009 was 5,556 Tg, with 89% of the C stored in soils. The estimated loss of C as CO2 and biogenic methane (CH4 ) emissions resulted in wetlands of Alaska increasing the greenhouse gas forcing of climate warming. Simulations for the projection period were conducted for six climate change scenarios constructed from two climate models forced under three CO2 emission scenarios. Ecosystem C storage averaged among climate scenarios increased 3.94 Tg C/yr by 2099, with variability among the simulations ranging from 2.02 to 4.42 Tg C/yr. These increases were driven primarily by increases in net primary production (NPP) that were greater than losses from increased decomposition and fire. The NPP increase was driven by CO2 fertilization (~5% per 100 parts per million by volume increase) and by increases in air temperature (~1% per °C increase). Increases in air temperature were estimated to be the primary cause for a projected 47.7% mean increase in biogenic CH4 emissions among the simulations (~15% per °C increase). Ecosystem CO2 sequestration offset the increase in CH4 emissions during the 21st century to decrease the greenhouse gas forcing of climate warming. However, beyond 2100, we expect that this forcing will ultimately increase as wetland ecosystems transition from being a sink to a source of atmospheric CO2 because of (1) decreasing sensitivity of NPP to increasing atmospheric CO2 , (2) increasing availability of soil C for decomposition as permafrost thaws, and (3) continued positive sensitivity of biogenic CH4 emissions to increases in soil temperature.

Fuel-reduction management alters plant composition, carbon and nitrogen pools, and soil thaw in Alaskan boreal forest.

Increasing wildfire activity in Alaska's boreal forests has led to greater fuel-reduction management. Management has been implemented to reduce wildfire spread, but the ecological impacts of these practices are poorly known. We quantified the effects of hand-thinning and shearblading on above- and belowground stand characteristics, plant species composition, carbon (C) and nitrogen (N) pools, and soil thaw across 19 sites dominated by black spruce (Picea mariana) in interior Alaska treated 2-12 years prior to sampling. The density of deciduous tree seedlings was significantly higher in shearbladed areas compared to unmanaged forest (6.4 vs. 0.1 stems/m2 ), and unmanaged stands exhibited the highest mean density of conifer seedlings and layers (1.4 stems/m2 ). Understory plant community composition was most similar between unmanaged and thinned stands. Shearblading resulted in a near complete loss of aboveground tree biomass C pools while thinning approximately halved the C pool size (1.2 kg C/m2 compared to 3.1 kg C/m2 in unmanaged forest). Significantly smaller soil organic layer (SOL) C and N pools were observed in shearbladed stands (3.2 kg C/m2 and 116.8 g N/m2 ) relative to thinned (6.0 kg C/m2 and 192.2 g N/m2 ) and unmanaged (5.9 kg C/m2 and 178.7 g N/m2 ) stands. No difference in C and N pool sizes in the uppermost 10 cm of mineral soil was observed among stand types. Total C stocks for measured pools was 2.6 kg C/m2 smaller in thinned stands and 5.8 kg C/m2 smaller in shearbladed stands when compared to unmanaged forest. Soil thaw depth averaged 13 cm deeper in thinned areas and 46 cm deeper in shearbladed areas relative to adjacent unmanaged stands, although variability was high across sites. Deeper soil thaw was linked to shallower SOL depth for unmanaged stands and both management types, however for any given SOL depth, thaw tended to be deeper in shearbladed areas compared to unmanaged forest. These findings indicate that fuel-reduction management alters plant community composition, C and N pools, and soil thaw depth, with consequences for ecosystem structure and function beyond those intended for fire management.

The role of driving factors in historical and projected carbon dynamics of upland ecosystems in Alaska.

It is important to understand how upland ecosystems of Alaska, which are estimated to occupy 84% of the state (i.e., 1,237,774 km2 ), are influencing and will influence state-wide carbon (C) dynamics in the face of ongoing climate change. We coupled fire disturbance and biogeochemical models to assess the relative effects of changing atmospheric carbon dioxide (CO2 ), climate, logging and fire regimes on the historical and future C balance of upland ecosystems for the four main Landscape Conservation Cooperatives (LCCs) of Alaska. At the end of the historical period (1950-2009) of our analysis, we estimate that upland ecosystems of Alaska store ~50 Pg C (with ~90% of the C in soils), and gained 3.26 Tg C/yr. Three of the LCCs had gains in total ecosystem C storage, while the Northwest Boreal LCC lost C (-6.01 Tg C/yr) because of increases in fire activity. Carbon exports from logging affected only the North Pacific LCC and represented less than 1% of the state's net primary production (NPP). The analysis for the future time period (2010-2099) consisted of six simulations driven by climate outputs from two climate models for three emission scenarios. Across the climate scenarios, total ecosystem C storage increased between 19.5 and 66.3 Tg C/yr, which represents 3.4% to 11.7% increase in Alaska upland's storage. We conducted additional simulations to attribute these responses to environmental changes. This analysis showed that atmospheric CO2 fertilization was the main driver of ecosystem C balance. By comparing future simulations with constant and with increasing atmospheric CO2 , we estimated that the sensitivity of NPP was 4.8% per 100 ppmv, but NPP becomes less sensitive to CO2 increase throughout the 21st century. Overall, our analyses suggest that the decreasing CO2 sensitivity of NPP and the increasing sensitivity of heterotrophic respiration to air temperature, in addition to the increase in C loss from wildfires weakens the C sink from upland ecosystems of Alaska and will ultimately lead to a source of CO2 to the atmosphere beyond 2100. Therefore, we conclude that the increasing regional C sink we estimate for the 21st century will most likely be transitional.

Can snowshoe hares control treeline expansions?

Treelines in Alaska are advancing in elevation and latitude because of climate warming, which is expanding the habitat available for boreal wildlife species, including snowshoe hares (Lepus americanus). Snowshoe hares are already present in tall shrub communities beyond treeline and are the main browser of white spruce (Picea glauca), the dominant tree species at treeline in Alaska. We investigated the processes involved in a "snowshoe hare filter" to white spruce establishment near treeline in Denali National Park, Alaska, USA. We modeled the pattern of spruce establishment from 1970 to 2009 and found that fewer spruce established during periods of high hare abundance. Multiple factors interact to influence browsing of spruce, including the hare cycle, snow depth and the characteristics of surrounding vegetation. Hares are abundant at treeline and may exclude spruce from otherwise optimal establishment sites, particularly floodplain locations with closed shrub canopies. The expansion of white spruce treeline in response to warming climate will be strongly modified by the spatial and temporal dynamics of the snowshoe hare filter.

Inland waters and their role in the carbon cycle of Alaska.

The magnitude of Alaska (AK) inland waters carbon (C) fluxes is likely to change in the future due to amplified climate warming impacts on the hydrology and biogeochemical processes in high latitude regions. Although current estimates of major aquatic C fluxes represent an essential baseline against which future change can be compared, a comprehensive assessment for AK has not yet been completed. To address this gap, we combined available data sets and applied consistent methodologies to estimate river lateral C export to the coast, river and lake carbon dioxide (CO2 ) and methane (CH4 ) emissions, and C burial in lakes for the six major hydrologic regions in the state. Estimated total aquatic C flux for AK was 41 Tg C/yr. Major components of this total flux, in Tg C/yr, were 18 for river lateral export, 17 for river CO2 emissions, and 8 for lake CO2 emissions. Lake C burial offset these fluxes by 2 Tg C/yr. River and lake CH4 emissions were 0.03 and 0.10 Tg C/yr, respectively. The Southeast and South central regions had the highest temperature, precipitation, terrestrial net primary productivity (NPP), and C yields (fluxes normalized to land area) were 77 and 42 g C·m-2 ·yr-1 , respectively. Lake CO2 emissions represented over half of the total aquatic flux from the Southwest (37 g C·m-2 ·yr-1 ). The North Slope, Northwest, and Yukon regions had lesser yields (11, 15, and 17 g C·m2 ·yr-1 ), but these estimates may be the most vulnerable to future climate change, because of the heightened sensitivity of arctic and boreal ecosystems to intensified warming. Total aquatic C yield for AK was 27 g C·m-2 ·yr-1 , which represented 16% of the estimated terrestrial NPP. Freshwater ecosystems represent a significant conduit for C loss, and a more comprehensive view of land-water-atmosphere interactions is necessary to predict future climate change impacts on the Alaskan ecosystem C balance.

Historical and projected trends in landscape drivers affecting carbon dynamics in Alaska.

Modern climate change in Alaska has resulted in widespread thawing of permafrost, increased fire activity, and extensive changes in vegetation characteristics that have significant consequences for socioecological systems. Despite observations of the heightened sensitivity of these systems to change, there has not been a comprehensive assessment of factors that drive ecosystem changes throughout Alaska. Here we present research that improves our understanding of the main drivers of the spatiotemporal patterns of carbon dynamics using in situ observations, remote sensing data, and an array of modeling techniques. In the last 60 yr, Alaska has seen a large increase in mean annual air temperature (1.7°C), with the greatest warming occurring over winter and spring. Warming trends are projected to continue throughout the 21st century and will likely result in landscape-level changes to ecosystem structure and function. Wetlands, mainly bogs and fens, which are currently estimated to cover 12.5% of the landscape, strongly influence exchange of methane between Alaska's ecosystems and the atmosphere and are expected to be affected by thawing permafrost and shifts in hydrology. Simulations suggest the current proportion of near-surface (within 1 m) and deep (within 5 m) permafrost extent will be reduced by 9-74% and 33-55% by the end of the 21st century, respectively. Since 2000, an average of 678 595 ha/yr was burned, more than twice the annual average during 1950-1999. The largest increase in fire activity is projected for the boreal forest, which could result in a reduction in late-successional spruce forest (8-44%) and an increase in early-successional deciduous forest (25-113%) that would mediate future fire activity and weaken permafrost stability in the region. Climate warming will also affect vegetation communities across arctic regions, where the coverage of deciduous forest could increase (223-620%), shrub tundra may increase (4-21%), and graminoid tundra might decrease (10-24%). This study sheds light on the sensitivity of Alaska's ecosystems to change that has the potential to significantly affect local and regional carbon balance, but more research is needed to improve estimates of land-surface and subsurface properties, and to better account for ecosystem dynamics affected by a myriad of biophysical factors and interactions.

Thermokarst rates intensify due to climate change and forest fragmentation in an Alaskan boreal forest lowland.

Lowland boreal forest ecosystems in Alaska are dominated by wetlands comprised of a complex mosaic of fens, collapse-scar bogs, low shrub/scrub, and forests growing on elevated ice-rich permafrost soils. Thermokarst has affected the lowlands of the Tanana Flats in central Alaska for centuries, as thawing permafrost collapses forests that transition to wetlands. Located within the discontinuous permafrost zone, this region has significantly warmed over the past half-century, and much of these carbon-rich permafrost soils are now within ~0.5 °C of thawing. Increased permafrost thaw in lowland boreal forests in response to warming may have consequences for the climate system. This study evaluates the trajectories and potential drivers of 60 years of forest change in a landscape subjected to permafrost thaw in unburned dominant forest types (paper birch and black spruce) associated with location on elevated permafrost plateau and across multiple time periods (1949, 1978, 1986, 1998, and 2009) using historical and contemporary aerial and satellite images for change detection. We developed (i) a deterministic statistical model to evaluate the potential climatic controls on forest change using gradient boosting and regression tree analysis, and (ii) a 30 × 30 m land cover map of the Tanana Flats to estimate the potential landscape-level losses of forest area due to thermokarst from 1949 to 2009. Over the 60-year period, we observed a nonlinear loss of birch forests and a relatively continuous gain of spruce forest associated with thermokarst and forest succession, while gradient boosting/regression tree models identify precipitation and forest fragmentation as the primary factors controlling birch and spruce forest change, respectively. Between 1950 and 2009, landscape-level analysis estimates a transition of ~15 km² or ~7% of birch forests to wetlands, where the greatest change followed warm periods. This work highlights that the vulnerability and resilience of lowland ice-rich permafrost ecosystems to climate changes depend on forest type.

Polygonal tundra geomorphological change in response to warming alters future CO2 and CH4 flux on the Barrow Peninsula.

The landscape of the Barrow Peninsula in northern Alaska is thought to have formed over centuries to millennia, and is now dominated by ice-wedge polygonal tundra that spans drained thaw-lake basins and interstitial tundra. In nearby tundra regions, studies have identified a rapid increase in thermokarst formation (i.e., pits) over recent decades in response to climate warming, facilitating changes in polygonal tundra geomorphology. We assessed the future impact of 100 years of tundra geomorphic change on peak growing season carbon exchange in response to: (i) landscape succession associated with the thaw-lake cycle; and (ii) low, moderate, and extreme scenarios of thermokarst pit formation (10%, 30%, and 50%) reported for Alaskan arctic tundra sites. We developed a 30 × 30 m resolution tundra geomorphology map (overall accuracy:75%; Kappa:0.69) for our ~1800 km² study area composed of ten classes; drained slope, high center polygon, flat-center polygon, low center polygon, coalescent low center polygon, polygon trough, meadow, ponds, rivers, and lakes, to determine their spatial distribution across the Barrow Peninsula. Land-atmosphere CO2 and CH4 flux data were collected for the summers of 2006-2010 at eighty-two sites near Barrow, across the mapped classes. The developed geomorphic map was used for the regional assessment of carbon flux. Results indicate (i) at present during peak growing season on the Barrow Peninsula, CO2 uptake occurs at -902.3 10(6) gC-CO2 day(-1) (uncertainty using 95% CI is between -438.3 and -1366 10(6) gC-CO2 day(-1)) and CH4 flux at 28.9 10(6) gC-CH4 day(-1) (uncertainty using 95% CI is between 12.9 and 44.9 10(6) gC-CH4 day(-1)), (ii) one century of future landscape change associated with the thaw-lake cycle only slightly alter CO2 and CH4 exchange, while (iii) moderate increases in thermokarst pits would strengthen both CO2 uptake (-166.9 10(6) gC-CO2 day(-1)) and CH4 flux (2.8 10(6) gC-CH4 day(-1)) with geomorphic change from low to high center polygons, cumulatively resulting in an estimated negative feedback to warming during peak growing season.